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Abstract: 

Let Xi = {Xi(t), t e T} be i.i.d. copies of a centered Gaussian process 
X = {X(t), teT} with values in WL d defined on a separable metric space T. 
It is supposed that X is bounded. We consider the asymptotic behaviour of convex 
hulls 

W n = conv{X 1 (t),...,X„(t), teT} 
and show that with probability 1 

lim 1 W„ = W 

(in the sense of Hausdorff distance), where the limit shape W is defined by the covari- 
ance structure of X: W — conv{Kt, t e T}, Kt being the concentration ellipsoid of 
X{t). 

The asymptotic behavior of the mathematical expectations Ef(W n ), where / is 
an homogeneous functional is also studied . 
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1 Introduction 

Let T be a separable metric space. Let Xi = {Xi(t), t <G T} be i.i.d. copies of a 
centered Gaussian process X = {X(t), t £ T} with values in M. d . Assume that 
X has a.s. bounded paths and consider the convex hulls 

W n = amv{X 1 (t),...,X n (t), teT}. (1) 

We are studying the existence of a limit shape for the sequence {W„}. 

Our work is motivated by recent papers U [5] inspired by an interesting 
implication in ecological context in estimating the home range of a herd of 
animals with population size n. Mathematical results of these articles consist in 
exact computation of a mean perimeter L n and area A n of W n in the case when 
d = 2 and X is a standard Brownian motion on T = [0, 1]. It was shown that 

L n ~ 27i"v / 2 Inn, A n ~ 27rlnn, n — > oo. (2) 

The relation between L n and A n being the same as the relation between the 
perimeter and area of a circle of the radius v21nn, it seems credible to sup- 
pose that W n rounds up with the growth of n. Our aim is to show that this 
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phenomenon really occurs for all bounded Gaussian processes. Our main result 
(Theorem 1) establishes the existence with probability 1 of the limit 



lim 



1 



(3) 



V21nn 



W n = W 



(in the sense of Hausdorff distance) and gives the complete description of the 
limit set W which is natural to call limit shape for convex hulls W n . In particular 
case of standard Brownian motion on [0, 1] the set W coincides with the unit 
ball B d (0,l) of W l . 

An interesting consequence of ((3|) is that the rate of the growth of the convex 
hulls W n is the same for all bounded Gaussian processes. 

The proof for continuous Gaussian processes may be easily deduced from 
the known results concerning the asymptotic of Gaussian samples (see [2 [5]), 
but in general case one needs an independent demonstration. 

Let us remark in addition that if T is a singleton, T = {to}) and d = 1, then 
the process X is simply a real random variable and W n is the segment 
[max{ X\, . . . , X n , } , min{ X\, . . . , X n , }]. It means that in some sense our study 
is closely connected with the classical theory of extrema. 

2 Asymptotic behavior of W n 
2.1 Notation 

Bd(0, 1), 5^(0, 1) are respectively unit ball and unit sphere of M. d . 
(•, •) is the scalar product in M. d . 

JC(B) is the space of compact convex subsets of a Banach space B provided with 
Hausdorff distance p^ : 



T is a separable metric space. 

C(T) is the space of continuous functions on T with uniform norm. 

X = {X(t), t € T} is a separable bounded centered Gaussian process with 

values in M. d . 

R t is the covariance matrix of X(t). 

K t is the ellipsoid of concentration of X(t) : 



pn(A, B) = max{inf{ e | A C B e }, inf{ e\B C A e }}, 



A e is the open e-neighbourhood of A. 
We set JC d = /C(R d ) and p = p Rd . 

M.a{9), S <G Sd(0, 1), is a support function of a set A E JC 



Ma{0) = sup (x, 9), 0eS d (O,l). 



K t = {x G R d | {R^ x x,x) < 1}. 



Finally we set 



W = conv{A" t , t G T}. 



(4) 
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2.2 Limit shape 
Theorem 1. 

1) Let X = {X{t), t G T} be a bounded centered Gaussian process with 
values in M d . Let (Xi) be a sequence of i.i.d. copies of X and W n be the convex 
hull defined by (Qp. 

Then with probability 1 

W n ^ W. (5) 



V21nn 

2) If T is compact and X is continuous, then a.s 



W n , w)=o(-±=). 



V2 hi n J \ \/m 

Remark 1. It is not difficult to see that the support function Aiw of the limit 
shape W admits the following representation 

M w {9) = cr{6), 

where 

cr 2 (6) =sup{R t 0, 9), 9eS d (0,l). 
teT 

The examples below show that in concrete cases the identification of W is not 
very complicated. 

Remark 2. For non-centered processes the relation ([5]) remains the same whe- 
reas ^ must be replaced by 

2.3 Asymptotic behavior of moments 

Let / : JC d — > R 1 be a continuous positive increasing homogeneous function of 
degree p , that is 

f{A) >0 Vie JC d - 

f(A 1 )<f(A 2 ) VA 1 cA 2 , A u A 2 €K d ; 
f(cA) = cPf(A), V c > 0, VA e JC d . 

Theorem 2. Let f be a function with the properties described above. Then, 
under hypothesis of Theorem 1 
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Remark 3. 

1) This theorem gives in particular the asymptotic behavior for mean values 
of all reasonable geometrical characteristics of W n (such as volume, surface 
measure, diameter, ...). 

2) By replacing / with f m , m > 0, we get the asymptotic behavior of higher 
order moments 

Ef m (~^= W n ) -> f m (W) . (9) 
\v2mn / 

2.4 Examples 

Brownian motion. Let X be a standard d-dimensional Brownian motion on 
T = [0,1]. Then K t = VtB d (0,l) and the limit shape is W = B d (0, 1). In 
particular, Theorem 2 gives for d = 2 the relations ([2|). 

Self-similar processes. Let X = {X(t), t 6 M + } be a Gaussian centered 
self-similar process (SSP) with values in K d . It means that for some a > the 
processes 

{X(at), t e {a a X(t), t G R+} 

have the same law for any a > 0. 

If we suppose that X is bounded, then we can apply our Theorem 1 to the 

restriction of X on [0, 1]. As X(t) = t a X(l), we have K t = t a K x which gives 
W = Kl 

This conclusion is available in particular when X has the stationary incre- 
ments: indeed, in this case the process is continuous as for any 6 G 5^(0, 1) the 
process (X(t), 9) is a fractional Brownian motion (FBM). 

Fractional Brownian Bridge. Now let us suppose that the coordinates of 
the process Y{t) = {Y x (t), . . . , Y d (t)} are independent FBM's: 

EYi(t) = 0, r(t, s) := EY^Y^s) = ±(t 2H + s 2H - \t - s\ 2H ). 

Then the conditional process related to the condition Y(l) = 0, which can be 
called Fractional Brownian Bridge, coincides in distribution with the process 

X(t) = {X 1 (t),...,X d (t)}, X i (t) = Y i (t)-r(t,l)Y i (l), t£[0,l}. 

It is clear that K t = a(i)B d (0, 1), where a 2 (t) = t 2H - \{t 2H + 1 - |1 - t\ 2H ) 2 . 
The function a 2 reaches its maximum at t — ^ and c^ax = Wh ~ \ ■ Finally we 
see that W = o- max B d (0, 1). 

3 Proofs 

Proving Theorem 1. The theorem is a consequence of two following lemmas. 
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Lemma 1. Let Y be a r.v. such that for all / y < ^ 

Eexp{~/Y 2 } < oo. 

.Let (Yfe) fee a sequence of independent copies ofY. Then with probability 1 

limsup — max{Yi, . . . , Y n } < 1. 
n v 2 In n 

Lemma 2. Lei Y = swp T X(t), where X = (X(t), t G T) is a centered bounded 
Gaussian process with sup T VarX(t) = 1 and let (Y n ) be a sequence of indepen- 
dent copies of Y. Then with probability 1 

lira inf max{ Y\ , . . . , Y n } > 1. 

" V 2 In n 

Proof of Theorem 1, first part. Fix 9 e 5^(0, 1). Define 7r e : R d -> K 1 by 

7Tg(a;) = (6,x), x <E M. d , and set 

MW n ) = [m$\MW], 

where 

= min{ mi,e}, mi,e = inf (9,Xi(t)), 

i<n t£LT 

AfW = msx{M it e}, M ifi = sup(#, JSQ(t)). 

Since the paths of X are bounded, due to the well known result of Fernique - 
Marcus and Shepp [3j [9] we have that 

EcxpijM^ 2 } < oo 

for all 7 < 2^(6) ' wnere & 2 {9) = sup teT Var(#, -Xj(i)). 
Then by Lemma 1 with probability 1 

limsup —t=M^ <a(0). (10) 
n v 2 In n 

On the other hand, by Lemma 2 

liminf - F =^Af< i e) > a (6), a.s. (11) 
" V21nn 

Therefore 

lim-^=^=^), a.s. (12) 
™ V21nn 

By the same arguments 

Hm-=|=mi 9 ) = -^), a.s. (13) 
" V21nn 



It means that for any 9 £ 5^(0, 1) with probability 1 



( ^=^« 



Let {ei, i = 1, . . . , d} be a basis of R d . Consider the parallelepiped C n defined 
by the orthogonal projections of -^==W„ onto coordinate axes. The relation 
(fli)) implies that 

d 

C n — > JJ[-o"(ei), c(ei)], a.s. 



Since —r= — W n C C n , the sequence {—r= — W n \ is bounded, and hence rela- 

V2 In n L V2 Inn J 

tively compact, in JC d . 

Due to the natural isometry between (K, d , p) and C(Sd(0, 1)) it follows from 
this that the sequence {M. n (9), 9 £ Sd(0, 1)} of support functions of the sets 
{-^==Wn} is a.s. relatively compact in the space C(Sd(0, 1)). Let 6 be a 
countable dense subset of 5^(0,1). Using the relation (|12p we see that with 
probability 1 for all 9 £ 6 

M n {9) ->■ <x{6). 



Together with relative compactness this shows that almost surely the sequence 
{■M n {-)} has a unic limit point. Then the same is true for {-j====W n }, and 
Remark 1 concludes the proof of the first part . 



Proof of Theorem 1, second part. Now we can consider the processes 
X and Xk as random elements of the separable Banach space B = C(T). By 
Theorem 2.1. of [5] with probability 1 

p K (W n , V21nnW?) = o(l), (15) 

where W n = convn{Xi, . . . ,X n } and W is the ellipsoid of concentration of X. 
Let ip : B -> /C d be defined by 

y?(a;) = conv{x(i), i e T}. 

It is clear that <p(W n ) = W„, (fi(W) = W, and it is easy to check that the map 
ip is Lipschitzian: 

p(<p(x), <p(y)) < pm(x,y), x,y £ B. 
Therefore (J5]) follows directly from (TTS")) . □ 



Proofs of Lemmas 1-2. 

Proof of Lemma 1. Let s > 0. Setting 

Z„ = = max{Yi , . . . , Y n } 
V2mn 
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our assumption implies 

P {Z„ > 1 + s} < nP{Y > (1 + s)V2\nn} < 

ngexp{ 7 F 2 } = s 
~ exp{ 7 (l + s) 2 21nn} U > ' 

where (5 = 1- 2j(l + s) 2 > if 7 > wlTTP- 
This inequality shows that the series 

^P{Z m a > 1 + S} 

m 

is summable if a > 1/8. We apply the Borel-Cantelli lemma and since s is 
arbitrary, we find that 

lim sup Z m a < 1 . 

m 

As for k € [m a , (m + l) a ) 



ln(m + 1) 



lnm 

we get 



lim sup Z n < 1 . □ 



PROOF of Lemma 2. Let < s < 1 and let t e T be chosen so that cr 2 := 
VarX(io) > s. We use the same notation Z n for — max{Fi, . . . , Y n }. We 
have 

P{Z„ <s} = P |y < sV21nn|" = F(sV2\nn) n 



< exp{-n(l - F(sV21nn))}, 

where F is the distribution function of Y. 
Note that 



1 - F(x) = P{swpX(t) > x} > P{X{t a )>x} > C^-^xpl-^V 
t I 2cr J 



which shows that 

It means that the series P {Z n _ s} is summable, and by applying the Borel- 
Cantelli lemma we finish the proof. □ 



Proving Theorem 2. The proof is based on the following lemma complet- 
ing the information given by Lemma 1. 



7 



Lemma 3. Let Y be a r.v. such that for some 7 > 

EexpilY 2 } < 00. 
Let (Yi) be a sequence of independent copies ofY. Then for any k G N 

( \ \ k 
supE —= max{ Y\,..., Y n } I < 00. 
n \ V 2 In n J 

Proof of Theorem 2. Due to the continuity of / and the convergence §5$ 
the result will follow from the uniform integrability of the family |/ ( J^™ ^ j . 
Using the notation from the proof of the first part of Theorem 1, we set 

m n = min m^, M n = max M^, D n = max{— m n , M n }, 

i=l,...,d i=l,...,d 

L n = [— D n , D n ] d . 

Since W n C L n , we have 

Kvffls) 5 = (^)" /(| - M1 "»' 

hence it is sufficient to state that for all p > 

/ D \ p 

sup <oo. (16) 



11 V V 2 ln n , 

The latter relation follows directly from Lemma 1 . and the theorem is proved. □ 



Proof of Lemma 3. Let 

Z n = = max{ Yi, . . . , Y n }. 
y2 Inn 

Denote by F, F n the distribution functions of Y and Z n respectively. By Markov 
inequality and by the assumption of Lemma 



1 - F n (x) < nP{Y > xV2hm} < An 1 ^ 2 ^ . (17) 
Hence for a = -L= 

r-oo r-oo 

E(Z n ) k = I x^il- F n (x))dx <a k + kAJ a^V" 2 ^ 2 dx. (18) 
n 1-27 ^ 2 < exp {— 7a; 2 Inn}, 



As for x > —h= 

— V2~f 
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we find that 

r°° , 2 

limsup / x k ~ 1 n 1 ~ 2 ' yx ~ dx = 0, 

n J a 

and we get from (fT8|) 

limsup£ , (Z n ) fc < a k , 

n 

which completes the proof. □ 



4 Concluding remarks 

1. It is clear that the result of the second part of Theorem 1 is still available 
if the space M. d is replaced by a separable Banach space. 

On the contrary, the similar question about the first part of Theorem 1 and 
about Theorem 2 is more delicate: their proofs are essentially based on the 
compactness of bounded subsets of M. d and it is not clear how to handle this 
obstacle in the infinite-dimensional case. 



2. The second interesting question is about the character of approach of 
i=B to W. Is it true that b n p ( , , W] does converge in law to some 

s/2 In n n l V V2 In n ' j ° 

limit for some choice of normalizing constants b n 1 

The same question can be also asked for the processes 
{b n (M n {9) - M(9)), e S"^ 1 }, where M n {6), M{9) are respectively the 
support functions of ^= = and W. It seems that the recent paper j6| may be 
useful in this context. 

3. What can we say on the behaviour of W n in non-Gaussian case? It is 
more or less clear that the convergence a.s. must be replaced by the weak one 
and the normalizing constants will be transformed from logarithmic to power 
ones. Indeed, if X is a vector in M. m with regulary varying distribution and (Xi) 
is a sequence of i.i.d. copies of X, then it is well known that the point processes 



converge weekly to some Poisson point process IIq,. 
It follows immediately from this that 



► conv(n Q ). (19) 



• I/O 



This fact is still available for a much more general case when X is a random 
element of an abstract convex cone K (see [5]) provided X satisfies the condition 
of regular variation (condition (4.5) in [2]). Hence the convergence (fT9)) may be 
considered as a "regular varying" analog of the second part of Theorem 1. The 
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main difficulty now is how to check the condition of regular variation for concrete 
situations. 
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